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I Abstract. Interacting galaxies like the famous Antennae (NGC 4038/4039) or 

Stephan's Quintet (HCG 92) show considerable star forming activity in their tidal 
(N ■ arms. High resolution images (e.g. from HST-observations) indicate that these re- 

^ , gions consist of up to hundreds of massive stellar clusters or tidal dwarf galaxies 

■ (TDG). In this paper we want to investigate the future fate of these clusters of 

— ' massive star clusters (in this work called super- clusters). We simulate compact super- 

, clusters in the tidal field of a host-galaxy and investigate the influence of orbital and 

' internal parameters on the rate and timescale of the merging process. We show that 

it is possible that such configurations merge and build a dwarf galaxy, which could be 
. an important mechanism of how long-lived dwarf satellite galaxies form. A detailed 

study of the merger object will appear in a follow-up paper. 



Keywords: methods: numerical - galaxies: interaction - galaxies: dwarf galaxies 
^ ■ galaxies: star clusters - star clusters: merging 



1. Introduction 



\^ • High resolution images from HST-observations of the Antennae galaxies 



(Whitmore et al. 1999, Zhang &: Fall 1999) show that the star forming 
regions there consist of up to hundreds of young (ages 3-7 Myr) , com- 
pact massive star clusters with dimensions of a few pc. The young star 
clusters in the Antennae have effective radii -Rgff = 4 pc and masses 
of the order 10^ - 10*^ Mq. These clust ers are not evenly distributed 
but themselves clustered in super-clusters spanning regions of several 
hundred pc in projected diameter and have concentrated cores, i.e. the 
cluster-density in the centre of the super-cluster is higher than in the 
outer parts. The richness of the super-clusters spans from groups of 
only a few to super-clusters containing hundreds of new star clusters. 
In other systems like Arp 245, Due et al. (2000) find a bound stellar 
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and gaseous object at the tip of the tidal tail. (Their numerical models 
show that the system is seen at approximately 100 Myr after the closest 
approach of the interacting pair.) This so-called tidal dwarf galaxy 
(TDG) contains old and new stellar material, a lot of gas and also 
undissolved massive stellar clusters. Systems like the merger remnant 
NGC 7252 (Miller et al. 1997) show a distribution of young massive star 
clusters which has the same age as the interaction [a 700 Myr). All 
these observations show that galaxy interactions lead to the formation 
of new massive and compact star clusters. 

Finally, galaxies like our Alilky Way host several dwarf galaxies 
(e.g. dSph or dE) (Mateo 1999) with high specific globular cluster 
frequencies (Grebel 2000) and globular clusters in their tidal streams. 

With this project we investigate the future fate of clusters of young 
massive star clusters. According to Kroupa (1998) it is possible that 
such configurations merge and build a dwarf galaxy. Therefore we sim- 
ulate compact super-clusters in the tidal field of a host-galaxy and 
investigate the influence of orbital and internal parameters on the rate 
and timescale of the merging process (i.e. how fast the single clusters 
merge and how many star clusters are able to survive this process) . In 
addition the properties of the resulting merger object and its dynamical 
evolution are studied. A detailed description of the properties of the 
merger objects will be given in a follow up paper. 



2. The Simulations 

The simulations are performed with the particle-mesh code SUPERBOX 
(Fellhauer et al. 2000). In SuPERBOX densities are derived on Cartesian 
grids using the nearest-grid-point scheme. Prom these density arrays the 
potential is calculated via a fast Fourier-transformation. The particles 
are then integrated forward in time using a fixed time-step Leap-Prog 
algorithm. SuPERBOX has a hierarchical grid architecture which in- 
cludes for each object two levels of high-resolution sub-grids. These 
sub-grids stay focused on the objects and travel with them through the 
simulation space, providing high resolution at the places of interest (in 
this case the super-cluster and the single clusters within). 

The massive star clusters are simulated as Plummer-spheres con- 
taining 100,000 particles each, having a Plummer-radius of Rpi = 6 pc 
and a cutoff radius i?cut = 30 pc. Each cluster has a total mass of 
= 10® M0 and a crossing time of 1.4 Myr. 

The super-cluster is also modelled as a Plummer distribution made 
up of A'^o star clusters described above, has a Plummer-radius and 
cutoff radius Rl^^^ = QR^i- The Plummer-radius of the super-cluster has 
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values of 50, 75, 150 and 300 pc. In this project the number of clusters 
is kept constant at Nq = 20, which is a typical number of star clusters 
found in these super-clusters. Nq was chosen small enough to get results 
with a considerable amount of CPU-time. Higher values of A^o will be 
dealt with using the newly available parallel version of SUPERBOX in 
the near future. In our calculations the super-clusters have an initial 
velocity-distribution according to the Plummer-distribution we gave 
them (i.e. they are initially in virial equilibrium). 

The super-cluster orbits through the external potential of a parent 
galaxy, which is given by 



with i?gai = 4 kpc and t;circ = 220 km/s. Wc refer the reader to Kroupa 
(1998), who describes the isolated case. The centre of the super-cluster 
moves on a circular orbit at distance D around the centre of the galaxy. 
The distance D from the galactic centre is varied to be 5, 10, 20, 30, 
50 and 100 kpc. 

The tidal radius Rt of the super-cluster depends mainly on D, but 
has also a low dependency on R^^- lies at the local maxima of $eff- 
It can be derived numerically by setting d^es/dr = where $efr is 



The grids in this project are chosen to have 64^ mesh-points with 
the following sizes: 

— The innermost grids cover single star clusters and have sizes (2 • 
Rcore) of 60 pc. This gives a resolution of 1 pc per cell. 

— The medium grids have sizes (2 • Rout) approximately equal to the 
cut-off radius of the super-cluster (i?cut) ensure that every star 
cluster is in the range of the medium grid of every other cluster. 
This means the medium grids have the sizes shown in Table I 

— The outermost grid (size: 2 • i?system) covers the orbit of the super- 
cluster around the galactic centre. This means 2 • -Rgystem is chosen 
to be 10 kpc larger than 2 • D, where D is the distance of the 
super-cluster from the galactic centre. 




(1) 
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Table I. Grid sizes and resolutions 
per cell of the medium grids. 



DSC 

-"•cut 


2 * Rout 


resolution 


300 pc 


600 pc 


10.0 pc 


450 pc 


1,000 pc 


16.7 pc 


900 pc 


2,000 pc 


33.3 pc 


1,800 pc 


3,000 pc 


50.0 pc 



The two-dimensional parameter-space of our simulations (R^i and 
D) can also be described with two dimensionless variables, namely 

a = Rpi/E^, (3 = RfjRt (3) 

a describes how densely the super-cluster is filled with star clusters. (5 
describes the strength of the tidal forces acting on the super-cluster. 

Table II lists the dimensionless parameters a and /3 for the different 
choices of the physical quantities i^^'j and D. 



Table II. Dimensionless parameters as function of the physical quantities. 



a 




psc 

^1 


[pc] 




P 




OSC 

-n-pl 


[pc] 




D [kpc] 


50 


75 


150 


300 


D [kpc] 


50 


75 


150 


300 


5 


0.12 


0.08 


0.04 




5 


0.775 


1.178 


2.500 




10 


0.12 


0.08 


0.04 




10 


0.617 


0.932 


1.935 




20 


0.12 


0.08 


0.04 


0.02 


20 


0.419 


0.630 


1.282 


2.761 


30 


0.12 


0.08 


0.04 


0.02 


30 


0.323 


0.486 


0.981 


2.048 


50 


0.12 


0.08 


0.04 


0.02 


50 


0.232 


0.348 


0.700 


1.431 


100 








0.02 


100 








0.89 



For each combination of (a, (5) several (3-6) random realisations 
are performed. Results discussed later for one pair of the parameter 
set are computed mean values out of the different simulations. Our 
study covers more than the observed ranges of parameters to show the 
influences of the different parameters more clearly. 

This is a first theoretical and numerical approach to investigate the 
fate of the super-clusters. Although, in nature, it is unlikely that two 
star clusters form overlapping but with Tincorrclatcd velocities, we do 
not reject random number placements which put two clusters in a 
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distance where they already overlap. Our simulations do not include 

gas dynamics, therefore they are valid after the expulsion of gas from 
the star clusters. By that time, clusters in the centre could already 
overlap each other. 



The number of star clusters in the super-cluster decreases due to two 
concurrent processes. The first and also the most important one is the 
merging process. The second one is the escape of star clusters from the 
super-cluster. Escape plays an important role only on long timescales 
or if the super-cluster is larger than its tidal radius (i.e. /3 > 1). 



where nm is the number of merged clusters and nesc the number of 
escaped clusters. 

3.1. Merging Timescales 

In our simulations the timescale of the merger process is very short. 
Most cluster merge within the first few crossing times of the super- 
cluster, forming a dense and spherical merger object in the centre of 
the super-cluster. 

To determine the timescale of this merger process we first have a 
look at simulations with (3 < 1.0. With this restriction we can neglect 
the number of escaping clusters (i.e. ncsc(i) = 0). 

Then wc take the following ansatz for the decrease of the number 
of clusters. If a star cluster makes one crossing through the super- 
cluster the chance to meet and merge with another cluster is the area 
covered by the cross-sections of all other clusters divided by the area 
of the super-cluster. We therefore define the merger-rate R (number 
of merger events per dimensionless time unit) as the ratio between 
the cross-sections of all (A^) star clusters (v4hit) travelling through the 
super-cluster and the area of the super-cluster A^c- The merger rate R 
per T^^ (crossing time of the super-cluster) should therefore be given 



3. Results 



N{t) = No - n^{t) - n^c{t) 



(4) 



by 



R = -'-^ 
dr 



^hit 



(5) 



with 



t 



(6) 



r = 



cr 
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being a dimensionless time and 

Ac = vr • (i?f,,,eff)' (7) 

is the projected area of the super-cluster with a mean radius that 
includes all clusters. This radius is smaller than i?^^^ due to the fact 
that the super-clusters contain only a limited number of clusters. In 
our simulations we find that 7 is given by7 = 3.8±0.1. 

For the cross-sections of all clusters we take the ansatz (see also 
Fig. 1) 

Ahit = (iV-l)-ylci with (8) 

^cl = TT • 




Figure 1. Schematic view of the first ansatz to determine the merging timescales. 

Every star cluster sees the cross-section A^i oi N — 1 other clusters 
and rm being the maximum distance at closest approach which leads 
to a merger afterwards. R is then proportional to N'^ and the number 
of clusters N(t) should decrease with time proportional to 1/(1 + kr). 

To determine ^ci we first have to check how the average encounter 
between two star clusters (or a star cluster and the merger object) looks 
like. Therefore we have to check the relative velocities and separations 
at the point of closest approach, which afterwards leads to the merging 
of the two star clusters. Fig. 2 shows the relative separation between a 
particular star cluster and the merger object. One can clearly see that 
the two objects are separated at their first encounter, implying the 
distance at closest approach is larger than the mean square radius of a 
single cluster. Their relative velocity is larger than the average velocity 
inside the super-cluster but less than a factor of three. Because even 
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Figure 2. A typical example for the merging of a star cluster with the merger 
object. Left panel: Solid line shows the distance between the merger object and 
the star cluster. Dashed lines show the half-mass- and the 90%-radii of the merger 
object. Horizontal line marks the maximum distance for the merger-criterion - if 
the distance of two objects stays less than this for the rest of the simulation the two 
objects are assumed to be merged; vertical line shows the adopted merger-time tm- 
Right panel: Ratio between the relative velocity Vrei of the cluster and the merger 
object and the velocity dispersion asc of the super-cluster. 



head-on encounters do not lead to merging if the encounter velocity lies 
above a certain treslihold, we follow a simple criterion extracted from 
Gerhard & Fall (1983 their Fig. 1) 

< — — (9) 
-Kpl 

where \/2crsc is the typical relative velocity of a pair of star clusters. 
This equation holds for a < 0.085. For a > 0.085 we argue below that 
in our parameter range we already start with a merged object in the 
centre of the super-cluster, and therefore M^i has to be replaced by 
uMci where u is the number of merged central clusters. With u > 2 
and increasing with increasing a, Eq. 9 holds for all values of a in our 
parameter space. 

In this paper, we define two clusters as merged if their mutual 
distance stays smaller then five Plummer-radii Rpi for the rest of the 
simulation. We found that 5 Rpi is a good compromise between declar- 
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ing no cluster as merged (the merger criterion being chosen too small) 
and merging all clusters right form the beginning (their mean distance 
being smaller than the merger criterion). Smaller merger radii have 
problems in the late stages of the simulation, where the centres of 
density of dissolved clusters are very hard to determine and could be 
found off-centre of the extended merger object. If the merger radius is 
chosen too small, these clusters would not be counted as merged. An 
energy criterion has to be handled with care, because clusters may be 
bound to the merger object without merging but staying on a circular 
orbit and decaying on the dynamical friction timescale. 

The next step is to approximate the energy gain of the clusters due 
to the passage. Aguilar & White (1985) showed that the total energy 
exchange in the tidal approximation given by Spitzer (1958), 

gives reasonable results if > br^- In this formula denotes the mean 
square radius of a single cluster. With a cut-off of 6i?pi we calculate the 
mean square radius of a Plummcr-sphere as approximately 1.5i?pi which 
2.3i?pi. Fig. 2 shows that the two objects pass each other for 
the first time at a distance which holds for the criterion of Aguilar k. 
White, r'p and Vp are the distance and the velocity at the point of closest 
approach. If we assume that Vp is equal to the mean relative velocity 
of two clusters in the super-cluster, y/losc, i-e. we neglect gravitational 
focusing and the acceleration of the star clusters firstly and taken into 
account that both clusters are able to gain energy we get 

4 • 2.3 G^M3i?2 
= 2^^ (11) 

To obtain the critical impact parameter which leads to a merger, we 
set the energy exchange equal to the orbital energy of the two clusters 
in the super-cluster: /S.E = j M^x (v^fsc)^- Inserting for the mean 
velocity of "particles" in a Plummer-sphere 

o 37rGM,e 



with Msc = NqMci and using the definition of a from Eq. 3, one obtains 
.322.2.3\'/' 



for Tp 



( 3(37r)2iV2 



Va-Rpi « 0.65Va-i?p1. (13) 
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In this formula we have not taken into account that the clusters are 
gravitationally focused. Inserting Spitzers (1987; eq. 6-15) formula for 
gravitational focusing, 



y rp(V2crsc)^ 



we obtain 



= 0.65 • ■ y 1.0 + ^ • R^^. (15) 

Given the values of a of our simulations we derive the values for rm 
given in Table III. The merger rate R is then 

Table III. Merger radius rm in units of iipi and -Rpi and 
proportionality factor 5No for the different a-values. 



a 


0.24 


0.12 


0.08 


0.04 


0.02 


0.006 






0.46 


0.36 


0.31 


0.25 


0.20 


0.14 






1.91 


2.97 


3.88 


6.18 


9.96 


23.4 


S{a 


)-iVo 


0.29 


0.18 


0.13 


0.087 


0.055 


0.027 



(0.65)2a(l + ^) 

R = N-(N-l)- ^ 6ia) ■ N^. (16) 

7 

Solving for R = — diV/dr gives 

Nir) = No ■ , (17) 

1 + o{a)f\oT 

with the values for (5(a) • Nq shown in Table III. This dependency is 
plotted as dashed lines in Figs. 3 & 7. 

On the other hand, using instead the merger criterion by Aarseth &; 
Fall (1980) as described in Gerhard & Fah (1983), 

2 9 

r V 
P I P < 1 ('-|o^ 

[4i?h]2 ^ [1.16^;esc(p)]2 - ' ^ ^ 

(their Eq. 4) where 

V iv? - (19) 



mfellhauer.tex; 5/02/2008; 5:40; p. 9 




Figure 3. Number of remaining objects vs. time measured in crossing times of the 
super-cluster for different values of the parameter a. (3 is kept smaller than 1.0. Lines 
show theoretical curves. Solid line corresponds to an exponential decrease, dashed 
line shows the 1/(1 + fcr)-decrease which depends on a. Dotted lines rely on the 
Aarseth & Fall criterion (Eq. 20). 

(their Eq. 5) with being the half-mass radius of a Plummer-sphere 
{k, 1.3i?pi) and = \/2crsc we obtain 

rp < (27.04 + 67.6A/671.75a2 + 22.31 • a + 1754.79 • a^) ^'"^ (20) 

Using Tp from this expression instead of in our Eq. 8 we obtain the 
dotted curves in Figs. 3 & 7. We note that this ansatz does not lead to a 
good description of our results. It leads to an a dependency of rp oc a^, 
which is steeper than we find in our simulations. Also the merger theory 
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of Makino & Hut (1997) does not agree with our results. They find a 
dependency on a with a power of 1.5, which does not fit our data. 
In our simulations we find a very weak dependency of the merging 
timescales on the parameter a. For our main parameter range (a = 
0.02 — 0.12) there is even no visible dependency. The simulations show a 
clear exponential decrease of the number of clusters with dimensionless 
time r independent of the choice of a (see Fig. 3). Therefore, a merger 
theory depending on the merging of single clusters must be wrong. 

The mean projected distance of clusters inside the innermost Plum- 
mcr radius of the super-cluster depends only on the choice of ATq and 
is given by 



<^mean(l-Rp'l) — <^ 



proj 



Numerical simulation shows that k, is given by 0.53. Since half the star 
clusters lie within one projected Plummer radius, we get for our choice 

of A?^3 = 20 ^ iVproj = 10 

Cimean(li?pl) ~ 0.307?^^. (22) 

Comparing this value with the merger radii rm for our choices of a (see 
Table III) one can see that in the centre of the super-cluster the mean 
distance between two clusters (within the innermost Plummer-radius 
i?p^) is comparable the merger radius for most a. In the centre the clus- 
ters are not able to separate from each other right from the beginning 
and the merging of them should therefore happen very quickly within 
one or two crossing times of the super-cluster, especially for high a. 

In this case we have to deal with a big merger object covering the 
central area of the super-cluster. Merging happens preferably with this 
central object. Therefore, our ansatz for A^it in Eq. 8 fails. Instead we 
use the following ansatz shown in Fig. 4. The merging cross-section 
Ahit is now a fraction e of the area of the whole super-cluster ^gc 

^hit = e^sc- (23) 

which has to be determined. Inserting in Eq. 5 (R = — dN/dr = eN) 
and integrating leads to 

N{t) = A^o •cxp(-er). (24) 

The exponential decrease, according to Eq. 24, is plotted in Fig. 3 
as the solid lines. There one sees clearly that the number of clusters 
first decreases exponentially. One should also notice that the merging 
timescale is now independent of the choice of a. 
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Figure 4- Schematic view of the new ansatz to determine the merging timescales. 

After all clusters travelling through the eentre have merged with the 
central merger object, the further decrease of N with time levels off. 
This is due to the fact that even if ^ < 1.0 there is the chance that 
some clusters escape and as a second effect there are clusters on rather 
circular orbits which do not travel through the central area. In low a 
cases (a = 0.04, 0.02) these clusters do not merge but their orbits will 
decay on timescales of the dynamical friction. 

To prove the independency on a we consider two characteristic 
quantities, namely the mean merger rate within the first crossing time 
(Ri.T^c) and the half- life merging time {T1/2), i.e. the dimensionless 
time until Nq/2 of the clusters have merged. Fig. 5 shows the number 
of clusters merging within the first crossing time, and the time it takes 
to merge half the clusters. All simulations with the same a, and which 
have /? < 1, are binned together. As can be seen, there is no dependency 
on a. 

The mean merger rate within the first crossing time and the mean 
half-life merging time derived from our simulations are given by 

i?i.r- = 3.90 ± 0.32 [merger events/T^^] (25) 
ri/2 = 4.18 ± 0.35 [T!^] 

are plotted in Fig. 5 and are also displayed in Fig. 6. Inserting this 
values in our ansatz (Eq. 23) for the merging times gives 

e 0.2. (26) 

Fig. 6 shows the mean values of Ri.t^c and T1/2 for all combinations 
of (a, P). It is clearly visible that there is no dependency of the results 
on /3 at all. It therefore seems that f3 influences the number of clusters 
which are able to merge rather than the timescale of the merger process. 
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Figure 5. Ri.t^^ and T1/2 as function of a for all simulations with /3 < 1. Solid lines 
show mean values from Eq. 25 while dotted line corresponds to the a-dependent 
theory from Eq. 16. 



We expect to obtain a change in the initial behaviour of the system 
if rjii falls well below the mean projected distance of the clusters. For 
our choice of A^O) this should happen approximately for a < 0.02. Then 
our first ansatz (Eq. 5) should represent the system. On the other 
hand if a becomes higher than ~ 0.2, the mean distance of all clusters 
is smaller than the merger radius - all A^o star clusters should then 
merge within one or two crossing times. To prove this behaviour we 
performed test calculations without tidal field. The results (Fig. 7) 
show nicely the transition between these three regimes. The upper left 
panel is a simulation with a = 0.24 and one sees clearly that almost 
all clusters merge within the first two crossing times. In the lower right 
panel a simulation with a = 0.006 is displayed. The number of clusters 
decreases according to our a-dependent theory. 



3.2. Merger-Rates 



The number of clusters, rim, which end up in the merger object shows 
no significant dependency on a and f3, as long as the super-cluster 
configuration is well inside it's tidal radius (i.e. P < 1.0). In this case 
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Figure 6. Ri.t^^ (boxes) and Tif2 (crosses) as a function of /3 for the different choices 
of a. SoUd hne is the mean value for T1/2, dashed line for Ri.t^^ from Eq. 25. 



almost all clusters merge and only one or two clusters sometimes survive 
by chance. As one can see in Fig. 8, Tim is close to A'^o = 20. 

This changes if (3 becomes larger than 1. There is a significant drop 
in Tijn which also shows a weak a-dependency. 

We interpret this result as follows: We start with Eq. 4 and neglect 
the rare case of an escaping merged cluster here. If /5 is small (/5 < 1) 
tidal effects are not dominant, and the evolution of N(t) is determined 
by the merger processes alone as discussed before. If /? > 1 there is 
a trend that clusters can leave the super-cluster before participating 
in the merger events. The number of escaping clusters depends on 
how many clusters initially are outside the tidal radius and on the 
individual velocities of such clusters. With only 20 clusters initially for 
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Figure 7. Same as Fig 3 but without tidal field and more extended Q-range. 
Errorbars are omitted because only one simulation for each Q-value was performed. 



the entire super-cluster any statistics of the subset of escaping clusters 
is extremely poor. Their number strongly depends on the random num- 
bers used for the initialisation of the system. Keeping in mind this poor 
statistical weight of our data, we nevertheless identify two physically 
reasonable trends in Fig. 8: first, the larger /3, the smaller the number of 
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Figure 8. Number of merged clusters rim vs. parameter /3 for different values of a. 



merging clusters, and second, this trend is more pronounced for small 
a. This result is consistent with the picture that strong tidal fields lead 
to more escapers (i.e. less clusters available for merging); if, however, 
the individual clusters are relatively extended (larger a), the merging 
competes with the escape; clusters on orbits of potential escapers could 
be captured in the central regions by merging with a higher probability. 

Binning the simulations in 5 /?-bins, as shown in Fig. 9, shows that 
the dependency on (3 is linear. As best-fit for the data of our simulations 
we calculate 

= 20 - (2.5 ± 0.1) • (3. (27) 

To advance that kind of reasoning to an extreme, we could compare 
our results with the findings of Baumgardt (1998), who models the 
escaping stars from star clusters. He proposes a /3^/^-dependency. Our 
results are consistent with Baumgardt's results, in particular for a ^ 
in which case N{t) — > A'o in Eq. 17 (where Uesc was neglected), so that 
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Figure 9. Number of merged cluster as function of /3. Simulations are put together 
in 5 /3-bins. Straight line is the fitted linear decrease; curved line would be the 
theoretical ''^-dependency as stated in Baumgardt (1998). 

N{t) ^ TVo ~ f^csc(i) (Eq. 4). Note, however, that our particle number 
is very small compared to that work. 

3.3. Building Up the Merger Objects 

The formation scenario of the merger object depends on the chosen 
a-value. One finds that with high values of a, one has overlapping 
star clusters at the centre right from the beginning, and the simulation 
already starts with a merger object in the centre of the super-cluster. 
With decreasing values of a, the merger-tree starts with the merging 
of clusters at different positions in the super-cluster. Afterwards, these 
merger objects sink to the centre and merge together. The exact details 
of the merging depend however on the starting conditions, and there 
can always be cases with high a behaving like low a and vice versa. 

Fig. 10 and 11 show some snapshots of the evolution of clusters with 
high and low a- values. 
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Figure 10. Surface-density contour plots (resolution: (10 pc)^) of a simulation with 
Of = 0.12 and D = 10 kpc. The crossing time of the super-cluster is 7.4 Myr. High 
values of a correspond to compact super-clusters with short crossing times. 



Since the crossing-time of the super-cluster expressed in physical 
units depends on its size and therefore on a, the merger rates, if 
expressed in Myrs, also depend on a: Large a-values correspond to 
compact clusters which merge within a few tens of a Myr, while small 
a-values correspond to extended clusters in which it can take up to 1 
Gyr until all clusters are finally merged. 

After 20 Myr almost every cluster has already fallen into the merger 
object at the high a simulation (Fig. 10). At t = 100 Myr only the last 
single cluster can be seen (as a disturbance in the contours) merging 
with the main object. 

In the low a simulation (Fig. 11) almost every cluster is seen as an 
individual object, even after 100 Myr. After 1 Gyr the merger object 
is surrounded by clusters which are still in the process of merging. 
These clusters may account for a high specific cluster frequency of 
the merger object. Not shown in the last snapshot are the 2 escaping 
clusters travelling on the same orbit around the galaxy as the merger 
object. 

The merger objects show an exponential density distribution with 
exponential scale-length rexp ~ 10 pc and follow a de-Vaucouleur sur- 
face density profile with line-of-sight velocity dispersions about 20 km/s. 
A detailed description and analysis of our merger-objects is subject of 
Fellhauer (2000) and a follow-up paper. Here we merely briefly note 
that our merger objects fall in the region between dwarf galaxies and 
globular clusters in the central surface brightness - absolute magnitude 
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Figure 11. Surface-density contour plots (resolution (10 pc)'^) of a simulation with 
a = 0.02 and D — 100 kpc. The crossing time of the super-cluster is 108.4 Myr. 



diagram (Fig. 3 in Ferguson & Binggeli 1994). Such objects evolve into 
dSph-like systems in a periodic tidal field (Fig. 13 in Kroupa 1997; 
Fig. 2 in Kroupa 1998b). More massive and extended merger-objects 
can, of course, be obtained by suitable choices of R^^ and Mgc, but we 
defer a more detailed discussion of this to the future. 



4. Conclusions 

We have performed a set of self-consistent dynamical models of clusters 
of twenty gas-free star clusters, as they have been recently observed 
in the Antennae galaxies. For all our models, a central merger object 
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formed out of a large fraction of the single clusters. As long as the 
super-cluster is smaller than its tidal radius, almost all clusters merge. 

For our choice of parameters (5 kpc < D < 100 kpc and 300 pc 
^ -^cut ^ 1-^ ^pc) the merger process does not depend significantly 
on the two main dimensionless parameters of the problem, the tidal 
field strength relative to the super-cluster concentration (/?) and the 
relative concentration of the super-clusters and the individual clusters 
(q). The cluster merging process cannot be modelled by a sequence of 
two-body merger events, but is a "collective" interaction, where the first 
passage ("encounter") of a cluster through the dense central region of 
the super-cluster leads to the assimilation into a growing merger object. 
The timescale of the merging process is the same (measured in internal 
crossing times of the super-cluster T^^) for all models. Measuring time 
in years shows that high a calculations form the merger object within 
very short times (~ 50 Myr), while for extended clusters (low a) this 
process can take up to 1 Gyr or even longer. 

While for very strong tidal fields the tidal mass loss as e.g. discussed 
by Baumgardt (1998) dominates and merging processes are suppressed, 
and for the limit where the individual clusters approach point masses 
the merging is suppressed as well, our parameter range, which is consis- 
tent with the observations, always allows for the quick merger scenario 
on a few crossing time scales. Tidal mass loss is only a secondary effect 
for part of our parameter space (stronger tidal field) due to a slight 
reduction of the number of clusters available for merging. The number 
of merged clusters nm decreases linearly with /3 if /? becomes larger 
than 1.0. 

We do not claim that our process of forming dwarf galaxies is the 
only possible way how these objects form, but at least it is one possible 
way to explain the existence of low-mass dwarf galaxies in the vicinity 
of large "normal" galaxies like our Milky Way. In addition, some "side- 
effects" can be explained by our models. It takes very long until all 
clusters finally end up in the merger object for low values of a. Such 
systems could account for a high specific globular cluster frequency 
found in dwarf spheroidal galaxies. Our models develop tidal tails which 
spread along the orbit of the dwarf galaxy. Surviving escaped star 
clusters are also found on the orbit. This is similar to what is observed 
for the Sagittarius dwarf spheroidal galaxy. 



Acknowledgements 

Part of this project (parallelisation of the code, use of high-performance 
computing power) was carried out at the Edinburgh Parallel Comput- 



mfellhauer.tex; 5/02/2008; 5:40; p. 20 



Merging Star Clusters 



21 



ing Centre (EPCC) by the TRACS-programme. The TRACS-program- 
mc is a scheme of the European Community (Access to Research In- 
frastructure action of the Improving Human Potential Programme; 
contract No HPRI-1999-CT-00026) which aUows young Phd- or post- 
doc scientist to learn parallel computing, offering computing time and 
support for their projects. MF thanks the staff of EPCC for their sup- 
port. We also acknowledge the insightful discussions with Prof. D.C. 
Heggie of the Department of Mathematics &; Statistics at the University 
of Edinburgh. 

References 

Aarseth, S.J., Fall, S.M., (1980) ApJ, 236 43. 
Aguilar, L.A., White, S.D.M., (1985) ApJ, 295, 374. 
Baumgardt, H., (1998) A&A, 330, 480. 

Binney, J., Tremaine, S., (1987) "Galactic Dynamics", Princeton Univ. Press. 
Due, P.-A., Brinks, E., Springel, V., Pichardo, B., Weilbacher, P., Mirabel, I.F., 

(2000) AJ, 120, 1238. 
Fellhauer, M., Kroupa, P., Baumgardt, H., Bien, R., Boily, CM., Spurzem, R., 

Wassmer, N., (2000) NewA, 5, 305. 
Fellhauer, M., Kroupa, P., (2000) ASP Conf. Ser., 211, 241.^ 

Fellhauer, M., (2001) PhD-thesis Ruprecht-Karls University, Heidelberg, Germany; 

Shaker, Aachen, Germany. 
Ferguson, H.C., Binggeli, B., (1994) A&ARev., 6, 67. 

Gallagher, S.C., Hunsberger, S.D., Charlton, J.C., Zaritsky, D., (2000) ASP Conf. 

Ser., 211, 247.^ 
Gerhard, O.E., Fall, S.M., (1983) MNRAS, 203, 1253. 
Grebel, E.K., (2000) ASP Conf. Ser., 211, 262.^ 
Hirashita, H., Kamaya, H., Takeuchi, T.T., (2000) PASJ, 51, 12. 
Hunsberger, S.D., (1997) BAAS, 29, 1406. 

Hunsberger, S.D., Gallagher, S.D., Charlton, J.C., Zaritsky, D., (2000) ASP Conf. 

Ser., 211, 254.^ 
Klessen, R.S., Kroupa, P., (1998) ApJ, 498, 143. 
Kleyna, J., Geller, M., Kenyon, S., Kurtz, M., (1999) AJ, 117, 1275. 
Kroupa, P., (1997) NcwA, 2, 139. 
Kroupa, P., (1998) MNRAS, 300, 200. 

Kroupa, P., (1998b) in "The Magellanic Clouds and other Dwarf Galaxies" eds. T. 

Richtler, J.M. Braun, Shaker, Aachen, Germany, (astro-ph/9804255). 
Makino, J., Hut, P., (1997) ApJ, 481, 83. 

Martinoz-Delgado, D., Apaicio, A., Gomez-Flechoso, M.A., (2000) to appear in 
Proceeding of the Euro conference "The evolution of Galaxies. I. Observational 
Clues", Granada May 2000. astro-ph/0009066 

Mateo, M., (1998) ARAA, 36, 435. 

Miller, B.W., Whitmore, B.C., Schweizer, F., Fall, S.M., (1997) AJ, 114, 2381. 

Spitzer, L., Jr., (1958) ApJ, 127, 17. 

Spitzer, L., Jr., (1987) "Dynamical Evolution of Globular Clusters", Princeton 

University Press. 
Struck, C, Phys. Rep., (1999) 321, 1. 



mfellhauer.tex; 5/02/2008; 5:40; p. 21 



22 



Fellhauer et al. 



Whitmore, B.C., Schweizer, F., (1995) AJ, 109, 960. 
Whitmorc, B.C., Schweizer, F., (1995) A,], 109, 1412. 

Whitmore, B.C., Zhang, Q., Leitherer, C, Fall, S.M., (1999) AJ, 118, 1551. 
Zhang, Q., Fall, S.M., (1999) ApJL, 527, 81L. 

^: ASP Conf. Ser. 211: Proceedings of the workshop "Massive Stellar 
Clusters" held in Strasbourg Nov. 1999; eds. A. Lancon, CM. Boily 



mfellhauer.tex; 5/02/2008; 5:40; p. 22 



